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ANALYSIS  OF  THE  DISPLACEMENT  FIELD  FOR  LOCALIZED  DISTURBANCES 
IN  A STRATIFIED  FLUID  WITH  SHEAR 


INTRODUCTION 

Many  theoretical  and  experimental  studies  have  been  devoted  to  various  aspects  of 
the  two-dimensional  collapse  of  a mixed  region  in  a stratified  fiuid.  Bell  and  Dugan  [1] 
and  Merritt  [2]  describe  work  in  this  area  up  to  about  1974.  Interest  in  the  phenomenon 
is  due  in  large  part  to  it*  relevance  as  a model  for  collapse  of  the  turbulent  wake  of  a self- 
propelhd  body  moving  through  a stratified  ocean. 

The  complexity  of  the  phenomenon,  which  involves  localized  turbulent  mixing  in  a 
stratified  fluid  and  the  generation  of  internal  waves,  has  led  many  investigators  (3-10]  to 
study  an  initial-value  problem  involving  a two-dimensional  disturbance.  Implicit  in  this 
simplification  is  the  assumption  that  no  mixing  occurs  following  the  initial  time.  Schooley 
and  Stewart  [11],  in  the  first  demonstration  of  the  wake  collapse  phenomenon  behind  a 
moving  self-propelled  body,  int-oduced  the  use  of  such  initial-value  problems  in  order  to 
gain  some  understanding  of  the  salient  features  of  the  collapse  process. 

We  present  here  an  analysis  of  an  initial-value  problem  involving  a localized  distur- 
bance in  a stratified  fluid  which  has  a background  shear.  Although  shear  flows  are  gen- 
erally present  in  the  ocean,  only  the  recent  work  of  Hartman  [6]  incorporates  shear  in  a 
mathematical  model  of  collapse.  Our  objective  is  to  demonstrate  that  curtain  features  of 
the  displacement  field  of  fiuid  particles  ran  be  calculated  directly  from  the  initial  data  of 
the  problem.  The  calculation  is  made  with  the  Lagrangian  viewpoint  and  is  valid  for 
finite-amplitude  disturbances  in  an  incompressible,  viscous  strati  lied  fluid  which  is  supposed 
nondiffusivc. 

In  addition  we  consider  the  inverse  question  whereby  inferences  about  tire  nature  of 
the  initial  disturbance  are  drawn  from  data  describing  the  shape  of  the  initial  perturbed 
region.  The  results  obtained  may  prove  useful  to  experimentalists  as  a means  for  esti- 
mating the  degree  of  homogeneity  of  mixed  regions. 

The  present  results  may  be  viewed  as  a generalization  of  earlier  work  (12]  in  wb  :h 
the  collapse  occurred  in  fluid  which  was  at  rest  outside  the  localized  disturbance.  From 
a theoretical  point  of  view  both  works  give  examples  of  the  fact  that  the  Lagrangian  view 
can  be  used  to  obtain  exact  information  in  some  rather  complicated  situations. 

An  elementary  application  of  the  results  is  made  to  simplify  an  argument  of  Stock- 
hausen, Clark,  and  Kennedy  (13]  concerning  zero  net  convergence  above  the  region  of 
initial  disturbance. 

Manuscript  submitted  August  16,  1976. 
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THE  INITIAL-VALUE  PROBLEM 
Case  of  Laminar  Motion 


The  mathematical  model  we  uae  is  that  of  an  incompressible,  viscous,  stratified  fluid 
which  is  regarded  as  nondiffusive.  The  assumption  that  the  fluid  is  nendiffusive  means 
that  those  quantities  affecting  the  density,  such  as  temperature  and  salinity,  are  not  altered 
by  molecular  diffusion.  Consequently  the  density  of  each  fluid  particle  remains  unchanged 
during  its  motion.  The  model  is  widely  used  in  the  study  of  stratified  fluids  [14]. 


The  Eulerian  form  of  the  governing  equations  is 


V*u  = 0 , 


dp 


ii  = °- 


a) 

(2) 


du 

dt 


= -Vp  + pg  + V*T 


(3) 


for  all  times  t > 0.  The  fluid  velocity  is  designated  oy  u,  the  density  by  he  pressure 
by  p,  the  constant  gravita*  onal  field  by  g,  and  the  viscous  stress  i ensor  by  T.  Tlie  time 
derivative  in  Eqs.  (2)  and  (3)  is  the  material  derivative. 


For  simplicity  we  assume  an  infinite  fluid,  with  no  free  surface  and  no  bottom 
boundary.  Although  the  analysis  to  be  carried  out  can  be  applied  to  many  situations 
involving  a free  surface  or  bottom  boundary,  the  discussion  of  all  possible  cases  would  be 
unnecessarily  long  and  might  conceal  the  simplicity  of  the  actual  considerations  involved 
in  ihe  analysis. 


We  use  two  coordinate  frames  (x,  y,  z)  and  ( X , Y,  Z)  that  are  coincident  to  describe 
the  motion  of  the  fluid.  A time-dependent  background  si  ar  motion  described  by  the 
velocity  :f?eld  (l/(y,  f).  0,  W(y,  t))  is  permitted.  Superposed  on  this  background  motion 
we  suppose  there  to  be  an  initial  localized  disturbance  describ’d  by  a local  velocity  dis- 
turbance and  a local  deformation  of  the  isopychnals,  both  independent  of  z.  Figure  1 is 
a sketch  of  .e  situation.  We  make  the  basic  assumption  that  the  background  motion  is 
stable  with  respect  to  the  initial  disturbance. 


Our  interest  is  in  the  motion  of  fluid  particles  x = x(X,  f).  defined  by 


x = X + d(X,f), 


(4) 


and  our  particular  interest  is  in  the  displacement  of  fluid  particles  d(X,  t).  Since  particles 
in  any  z = constant  plane  experience  the  same  forces,  the  displacement  field  d is  independ- 
ent of  Z. 


We  assume  that  at  t = 0 the  localized  disturbance  can  be  expressed  as 

p(X.  Y,  0)  = pe(Y)  + 5 p(X,  Y) 
and  the  initial  velocity  field  as 


(5) 
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\ REGION  OF  INITIAL  PERTURBATION 
\ 

\ 

\ 

\ 

\ 


ISORVCHNALS 


Fit- 1— Infinite  fluid  with  timrr-dependent  background  shear  flow  and  an  initial 
localised  disturbance.  t/(y.  I)  and  l?(y,  t)  are  components  of  the  shear  flow. 


U(X,  Y,  0)  = U(V,  0)  + 5tl(X,  Y) . ( 

The  finite-amplitude  initial  disturbance  described  by  6p  and  5u  is  assumed  to  vanish  as 
I A' I -*•  °°.  The  density  represents  the  stratification  prevailing  at  large  horizontal 
distances  from  the  localized  disturbance. 

Fror*  the  Eulerian  equations  we  find  that  the  backgroiu  d motion  must  satisfy 

b2U 

p‘  ( 


bW  _ b2W 
Pe  bt  ■ P ay2 


where  p is  the  mechanical  viscosity  of  the  fluid.  Equations  (7)  and  (8)  state  that  the 
background  motion  is  modified  by  diffusion.  Since  we  shall  deal  with  an  infinite  fluid, 
we  regard  the  background  motion  as  a specified  solution  of  Eqr.  (7)  and  (8). 

This  completes  our  description  of  the  formal  initial-value  problem  for  ’s.mina c flows. 
Before  proceeding  with  its  study,  we  first  make  some  remarks  i.bout  a similar  initial-value 
problem  which  ar  ses  in  the  case  of  localized  regions  of  turbulent  mixing  produced  by  the 
motion  of  a self-propelled  body  through  a stratified  fluid. 
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Remarks  on  the  Disturbance  Created  by  the  Motion  of  a Self-Propelled  Body 

We  give  a brief  discussion  of  a mathematical  model  for  the  generation  of  internal 
waves  by  the  turbulent  wake  of  a self-propelled  body.  The  model  originated  with  the 
work  of  Schooley  and  Stewart  [11]  and  has  received  much  attention  from  theorists  and 
experimentalists  [1-13].  Our  aim  is  to  demonstrate  that  the  governing  equations  for  the 
mean  disturbance  (the  ensemble  average  disturbance)  bear  a clo~e  resemblance  to  those 
for  the  laminar  flow  already  discussed  and  that  useful  results  follow  from  this  fact. 

When  a self-propelled  body  moves  rapidly  through  a stratified  ocean  environment, 
its  turbulent  wake  mixes  the  fluid  and  initially  increases  in  cross  section  through  entrain- 
ment of  fluid  outside  the  wake.  Eventually  the  intensity  of  the  turbulence  diminishes  to 
such  an  extent  that  buoyancy  forces  the  mixed  fluid  back  toward  an  appropriate  equilib- 
rium level  and  the  vertical  dimension  of  the  wake  reduces  in  size.  This  phase  of  the  wake 
behavior  is  referred  to  as  the  gravitational-collapse  phase  and  is  accompanied  by  the  prop- 
agation of  internal  waves  array  from  the  wake  region.  A primary  assumption  of  the  present 
mathematical  model  is  that  significant  mixing  of  the  fluid  no  longer  occurs  once  the  gravi- 
tational collapse  has  begun. 

Since  the  rate  at  which  the  moving  body  generates  mixed  fluid  is  usually  large  in 
comparison  to  the  rate  at  which  gravitational  collapse  occurs,  the  assumption  is  made  in 
the  model  that  the  mean  or  ensemble  average  characteristics  of  the  collapse  process  can 
be  regarded  as  independent  of  distance  along  the  track  of  the  body. 

We  express  the  disturbance  as 


U(x,  y,  t)  + u' , 

(9) 

P„(x,  y,  0 + p , 

(10) 

P(x , y,  t)  + p , 

(11) 

where  pOT,  U,  P are  ensemble  average  quantities.  We  use  (->  to  indicate  an  ensemble  aver- 
age; consequently  the  fluctuations  satisfy 

(“’>  - 0 • 

(o')  - 0 • 

<p'>-  o. 

Time  t is  measured  from  the  initiation  of  tne  collapse  process. 

The  ensemble  average  quantities  contain  the  mean  background  shear  flow,  the  mean 
density  anomaly  due  to  prior  turbulent  mixing:  of  the  fluid,  and  the  mean  internal  wave 
disturbance  resulting  from  cjllapse  of  the  density  anomaly.  The  basic  equations  govern- 
ing the  evolution  of  the  ensemble  average  disturbance  are 
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7-U  = 0, 


(p  Jf)  - -VP  + />„«  ♦ V-T 


The  absence  from  the  right  side  of  Eq.  (13)  of  term*  involving  the  fluctuating  quantities 
reflects  our  assumption  of  Q muring  during  the  collapse  phase. 

We  can  interpret  Eqs.  (12)  through  (14)  as  describing  the  motion  of  parties  which 
move  about  in  accord  with  the  ensemble  average  motion  U(x,  y,  f )•  Parties  paths  for  the 
ensemble  average  motion  : e calculated  from 


= U. 


Equations  (12)  and  (13)  correspond  to  Eqs.  (1)  and  (2)  for  the  laminar  motion  discussed 
earlier. 

It  shall  he^caie  clear  that  our  results  for  net  displacements  in  the  laminar  case  can 
also  be  deduced  for  the  ensemble  average  motion  just  described.  The  reason  for  this  is 
that  the  results  do  not  depend  on  the  explicit  form  of  the  momentum  equation,  only  on 
its  qualitative  features.  The  features  required  are  that  localized  disturbances  eventually 
decay,  thus  allowing  the  disturbed  fluid  to  return  to  an  equilibrium  state  specified  by  the 
background  shear  and  density  field. 


TKE  DIRECT  PROBLEM 
Preliminaries 

We  now  return  to  the  initial-value  problem  for  the  laminar  flow  outlined  earlier.  In 
the  direct  problem  the  initial  configuration  is  assumed  known.  The  aim  here  is  to  demon- 
strate how  certain  features  of  the  motion  can  be  deduced  without  having  to  solve  for  the 
complex  flow  field.  The  Lagrangian  view  is  used. 

On  physical  grounds  it  is  clear  that  as  time  progresses,  the  energy  contained  in  the 
localized  disturbance  will  be  radiated  into  the  surrounding  fluid  and  ultimately  dissipated 
by  viscous  forces.  Thus  eventually  t’4e  .notion  will  again  be  horizontal  in  accord  with 
Eqs.  (7)  and  (8),  with  p,(y)  descrying  the  stratification  throughout  the  fluid.  At  suf- 
ficiently large  times  the  velocity  of  particles  can  be  written  in  the  form 

^ = U(y*,t),  (15) 


£ = 0 
dt  u’ 
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dx 

dt 


where 


W(y*.  t) , 

(17) 

y*(X,Y) 

(1C) 

specifies  the  final  vertical  level  of  the  fluid  particle  initially  cl  X,  Y. 


The  forms  of  Eqi.  (15)  through  (17)  iiigyst  that  we  write  the  particle  motion  in  the 
forms 


X - x * y,  f)  + 

| Uiy*,  t)dt , 

(19) 

y = y + iHX,  y,  t) , 

(20) 

a = z ♦ r<X.  y,  t)  ♦ 

W(y*,t)dt. 

(21) 

If  there  were  no  initial  disturbance  to  the  fluid,  tlve  particle  motion  would  be  specified 
by  Eqs.  (19)  through  (21)  with  the  £,  t),  and  ? terms  absent.  Consequently,  £,  n,  and  f 
specify  the  full  contribution  to  the  displacement  field  arising  as  a consequence  of  the 
gravitational  collapse  of  the  initial  localized  disturbance 

From  Eqs.  (15)  through  (17)  we  conclude  that  at  large  times  (large  in  comparison 
to  t'Jte  effective  collapse  time,  which  is  controlled  by  the  Brunt- Viiisilii  period'  we  must 
have* 


lim  {(X,  V.  t)  = **(X,  V) , (22) 

( -+oo 

lim  T)(X,  y,  t)  = q*(X.  y) , (23) 

{-too 

lim  r(X,  y,  t)  = r*(X,  y) . (24) 

t -400 

We  shall  demonstrate  in  the  following  that  it  is  possible  to  compute  (*  and  tj*  directly 
from  the  initial  data  of  the  problem. 

We  are  making  some  formal  assumptions  here.  We  assume  that  the  initial-value 
problem  possesses  c well-behaved  solution  for  0 < t < °°.  (A  rigorous  proof  does  not 
exist.  1 This  requires,  in  part,  that  the  background  flow  must  be  stable,  and,  from  stabil- 
ity theory,  this  places  requirements  on  the  Richardson  number,  as  discussed  by  Turner  (15, 
pp.  100  If]. 

We  must  also  have  the  background  flow  be  stable  with  respect  to  the  localized  dis- 
turbance of  finite  amplitude  which  we  prescribe  through  the  initial  conditions.  The  collapse 
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may,  in  some  cates,  produce  internal  waves  which  interact  with  the  background  shear  in 
such  a fashion  as  to  lead  to  amplification  of  the  waves  (critical  layer  behavior)  and  to  the 
ponbiliiy  of  wave  breakdown  (15,  p.  124] . Such  behavior  seems  to  imply  that  our 
initial-value  problem  does  not  always  potseas  a well-behaved  solution;  however  general 
criteria  do  not  exist  In  a recent  paper  Hartman  (16]  treats  a linearised  initial-value 
problem  for  a twodimensional,  unbounded,  exponentially  stratified,  plane  Couette  flow. 
No  difficulties  are  encountered  because  of  critical  layer  behavior. 


Net  Vertical  Displacement 

From  Eq.  (2)  the  density  of  a particle  at  time  t must  be  the  same  as  the  den  sit}  of 
the  particle  initially.  This  yields 

Pl*(0.  y(0?  z(0.  t ] = pr(F)  + 6p(X,  Y) . (25) 

But  as  t -*  «,  y(f)  -*  y*(X,  F),  and,  since  the  ultimate  density  field  b again  described  by 
p*,  Eq.  (25)  becomes 


Pe{y*)  = Pe(Y  + rim)  = p,(F)  + 6p(X,Y).  (26) 

Equation  (26)  is  tlie  determining  condi tir*n  for  the  net  vertical  displacement  r}*{X,  V)  in 
terms  of  the  initial  data. 

Although  this  impucit  condition  for  tj*  can  be  solved  explicitly  for  many  stratifica- 
tions of  interest,  we  record  here  the  solution 

t,*(X  Y)  = (^)  6p(X,  Y) , (27) 

which  is  valid  for  those  cases  in  which  pe{Y)  is  linear  over  the  vertical  extent  of  the  initial 
disturbance. 


Net  Horizontal  Displacement 

For  the  determination  of  £*(X,  Y)  we  use 

d(x,  y.  Z) 

9(X  Y,  Z) 

which  is  a consequence  of  «n compressibility  and  conservation  of  mass.  The  notation 
indicates  the  Jacobi*:,  of  the  transformation  expressed  by  Eqs.  (19  through  (21).  Since 
the  x and  y coordinates  of  a particle  do  not  depend  on  Z , Eq.  (28  j can  be  simplified  to 


the  Lagrangian  condition 
= 1 , (28) 
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a* 

a* 

n 

ax 

ay 

u 

3x 

3x 

3y 

by 

A 

ax 

ay 

ax 

ay 

U 

by 

dz 

bz 

1 

ax 

ay| 

ax 

ay 

(29) 


For  the  class  of  problems  under  consideration  Eq.  (29)  places  no  restriction  on  ? and  our 
methods  will  yield  no  information  about  this  component  of  the  displacement  Held. 

Substitution  of  the  expressions  in  Eqs.  (19)  through  (21)  into  Eq.  (29)  yields,  when 
written  out. 


It  is  preferable  to  write  this  result  in  the  form 


ay  ax + 


U(y*,t)dt 


3y 

3X 


1, 


and  carrying  out  the  differentiation  of  the  integrals  provides 

'<  3 U{y*,t) 


an/,  a«  ^ r 
l1  + a x/v  * av/ " ay  ax  + [I 


dt 


( (30) 

\ax  ay  ay  ax/  1 • (d0) 


Because 


(by* 

dy* 

bjl  - 

\ ax 

ay  ' 

ay 

ax/: 

ax 

ay 

ay 

ax 

= 0, 


the  limiting  form  of  Eq.  (30)  as  t -*■  °°  becomes 

\ ay  / ax  ax  ay  ay ' 


(31) 


We  have  already  seen  how  rj*  can  be  deterniined  from  the  initial  data.  Equation  (31) 
then  represents  a linear  fust-order  partial-differential  equation  for  the  determination  of 
$*(X.  Y). 
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Comparison  to  the  Case  of  No  Shear 

The  equations  for  the  determination  of  p*(X,  V)  and  rj*(X,  V)  are  precisely  the  same 
equations  as  found  earlier  [12]  in  the  case  of  no  background  motion.  In  the  present  case 
£*  has  a different  interpretation,  since  the  shear  motion  itself  produces  a contribution  to 
the  displacement  field.  From  Ref.  12  we  have  at  once  that  £*  can  be  solved  for  by  using 
the  method  of  characteristics,  with  the  characteristic  equations  given  by 

dt*  dr?* 

*r  = -3y’  (32) 


do  37 


The  parameter  a identifies  points  on  a characteristic.  We  further  have  that  the  character- 
istics in  the  XY  plane  are  those  curves  specifying  the  initial  forms  of  ti’i  isopychnals. 

In  the  earlier  case  of  no  shear  [12]  we  confined  ourselves  to  symmetric  disturbances. 
We  then  had  £*(C,  Y)  = 0 as  initial  data  for  the  solution  of  Eqs.  (32)  through  (34).  In 
the  presence  of  shear  however  we  no  longer  have  such  a condition  even  for  symmetric 
disturbances.  Consequently  the  initial  data  available  allow  the  determination  of  only  the 
relative  quantity  p*(X,  Y)  - £q,  where  Pq  >s  the  unknown  valae  of  £*(A',  Y)  for  some 
point  on  the  same  isopychnal  that  the  point  (X,  Y)  is  on. 


A Simple  Application  to  the  Net  Convergence  of  Fluid  Particles 

Stockhausen,  Clark,  and  Kennedy  [13]  describe  a laboratory  experiment  in  which  a 
self-propelled  body  moved  through  a stratified  fluid.  They  give  an  argument  [13,  pp.  72 
and  78]  that  a net  convergence  between  two  fluid  particles  can  occur  only  when  the 
particles  are  at  a depth  level  where  mixing  has  taken  place. 

We  give  here  a rigorous  proof  that  the  statements  of  Stockhausen,  Clark,  and  Kennedy 
apply  to  the  initial-value  problem  we  have  been  examining.  Consider  two  fluid  particles 
P and  Q initially  on  an  undisturbed  isopychnal  as  shown  in  Fig.  2.  We  have  5 p(X,  7)  = 0 
for  all  points  on  the  isopychnal,  and  Eq.  (26)  provides  the  result  that  r?*(X,  7)  = 0 for 
all  points  on  the  isopychnal.  Using  this  result  in  Eqs.  (32)  through  (34),  we  get  cf£*/dX  = 0, 
aw*  this  can  bo  integrated  to  obtain  £*(P)  = £*(<3). 

Since  y*(X,  Y > = 7 for  all  points  on  the  isophchnal,  Eqs.  (19)  through  (21)  provide 


°,t)  = XP  + $(P,  t)  + [‘  U(Yp,  t)  < 

^ r\ 


and,  since  yQ  = YP, 
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Fig.  2— State  of  disturbance  at  t * 0 in  the  ay  plane,  with  the  dashed 
lines  indicating  the  region  containing  the  initial  disturbance.  Pand  Q 
identify  two  fluid  particles  tying  initially  on  an  undisturbed  iso- 
pychnal.  The  component  of  the  background  shear  flow  normal  to  the 
xy  plane  is  not  drawn. 


x(Q,  0 - XQ  + HQ,  t)  + 


Subtracting  the  two  expressions,  we  get 

x(P,  t)  - x (Q,  t)  = XP  - XQ  + HP  t)  - HQ,  t) . 

The  limit  of  this  as  f -*  °°  provides  the  ultimate  separation  of  the  particles  P and  Q in  the 
x direction  as 


x*(P)  - x*(Q)  = XP-  XQ  + ?(P)  - **(Q) 

= XP  - XQ, 

which  is  precisely  equal  to  their  initial  separation. 

We  conclude  that  no  net  convergence  occurs  between  particles  such  as  P and  Q lo- 
cated on  an  isopychnal  outside  the  region  of  initial  disturbance.  A net  convergence  can 
occur  only  if  P and  Q are  at  the  level  of  the  initial  disturbance  region.  This  is  perhaps 
an  unexpected  result.  One’s  initial  expectation  is  that  as  the  initial  disturbance  resolves 
itself  into  a system  of  internal  waves  which  propagate  through  the  body  of  the  fluid,  the 
end  result  would  be  particles  whose  locations  have  been  shifted  about  relative  to  one 
another.  The  severe  constraint  conferred  on  the  fluid  by  stratification  prevents  such  relative 
shifting. 
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REMARKS  ON  THE  INVERSE  PROBLEM 

A matter  of  great  practical  importance  is  the  nature  of  the  disturbance  in  the  fluid 
at  the  end  of  a brief  period  of  localized,  turbulent  mixing.  The  complexity  of  the  phe- 
nomena makes  such  information  difficult  to  obtain  by  direct  measurement. 

We  consider,  within  the  context  of  our  initial-value  problem,  an  indirect  method  for 
obtaining  information  about  the  initial  deformation  of  the  isopychnals.  In  effect  we  con- 
sider an  inverse  problem  in  which  it  is  assumed  that  the  background  stratification  is  known 
as  well  as  the  shapes  of  the  disturbed  region  at  t = 0 and  at  some  time  t large  in  compari- 
son to  the  time  scale  of  collapse.  An  earlier  report  [12]  dealt  with  this  question  in  the 
case  of  a zero  shear  environment.  Although  the  earlier  results  can  be  generalized  to  the 
situation  of  collapse  in  a shear  environment,  we  prefer  to  give  instead  a new  result  which 
seems  to  be  more  generally  useful. 

Consider  the  initial  configuration  of  the  fluid  as  shown  in  Fig.  3a.  By  hypothesis  we 
know  the  boundary  C within  which  the  initial  disturbance  is  confined.  A typical  iso- 
pychnal  of  density  pe  which  passes  through  C has  a form  as  shown  in  Fig.  3a.  Outside  C 
the  isrpychnal  is  horizontal  at  its  equilibrium  level,  while  inside  C it  is  deformed.  Let 
A(pe , t)  denote  the  area  enclosed  between  the  material  particles  initially  composing  the 
isopychnal  pe  and  the  particles  which  initially  form  the  boundary  C lying  above  the  iso- 
pychnal (Fig.  3).  As  a consequence  of  the  incompressibility  of  the  fluid  this  area  remains 
constant  throughout  the  motion,  so  that 


A(pe,t)  = A(pe,  0) 


for  all  time. 


We  now  prove  that  the  area  A(pe , t)  is  always  bounded  by  C and  the  isopych  pe. 
From  the  nondiffusive  nature  of  the  fluid  we  know  that  the  fluid  particles  retain  their 
density  throughout  the  motion.  Consequently  isopychnals  are  material  curves,  as  is  the 
curve  C.  We  conclude  that  A(pe,  t)  is  always  bounded  by  the  particles  comprising  C and 
the  isopychnal  pe. 

By  hypothesis  we  know  the  initial  and  “final”  configurations  of  the  material  contour 
C.  A final  configuration  of  the  fluid  is  one  in  which  the  isopychnals  must  be  horizontal 
and  at  the  same  level  as  given  by  the  known  initial  data.  Therefore  we  can  evaluate 
A(pe , °°)  for  any  isopychnal  passing  through  C and,  by  the  above  argument,  we  then  have 
A(pe,  C)  tor  each  isopychnal  passing  through  the  initial  contour  C. 

From  Fig.  3a  it  is  clear  that  knowledge  of  A(pe,  0)  can  be  used  to  calculate  the  area 
enclosed  between  the  position  of  the  deformed  isopychnal  pe  and  its  final  equilibrium 
level.  This  information  immediately  provides  the  average  level  of  the  isopychnal  pe  within 
the  initial  region. 

We  thus  have  a means  for  using  the  initial  and  final  shapes  of  the  disturbed  fluid  to 
estimate  average  characteristics  of  the  initial  disturbance.  This  method  may  prove  to  be 
of  use  in  estimating  such  quantities  as  the  degree  of  homogeneity  of  localized  disturbances. 
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(b)  / >>  time  *c«le  of  collapee 

Fig.  3— Initial  and  Anal  configurations  of  the  region  contain- 
ing the  initial  disturbance.  C denotes  the  boundary  of  the 
region  in  the  xy  plane.  A(pt,  t)  is  the  area  at  time  t enclosed 
between  the  isopychnal  of  density  pe  and  the  boundary  C. 


CONCLUSION 

We  considered  an  initial-value  problem  involving  a two-dimensional  localized  disturb- 
ance in  an  incompressible,  viscous,  stratified  fluid.  A time-dependent  background  shear 
flow  was  incorporated  in  the  analysis.  Although  the  problem  was  considered  within  the 
fra  u.ework  of  laminar  flow  theory,  we  pointed  out  that  similar  results  can  be  obtained 
for  the  ensemble-average  motion  resulting  from  the  collapse  of  a localized  region  of 
turbulent  mixed  fluid. 

The  results  are  perhaps  best  thought  of  as  an  extension  of  earlier  results  [1.2]  ob- 
tained for  the  case  of  no  background  motion.  Again  the  Lagrangian  view  has  been  used 
to  advantage. 

It  was  demonstrated  for  the  general  class  of  initial-value  problems  considered  that 
certain  components  of  the  net  displacement  field  could  be  calculated  directly  from  the 
initial  data.  This  calculation  is  effectively  exact  and  holds  for  finite-amplitude  disturb- 
ances. A simple  application  was  made  to  the  net  convergence  of  fluid  particles.  A poten- 
tial application  of  the  results  is  in  the  partial  verification  of  computer  algorithms  devised 
for  the  solution  of  initial-value  problems  of  the  type  considered. 
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We  also  examined  the  inverse  problem,  which  is  of  particular  interest  in  the  experi- 
mental investigation  of  localized  mixed  regions  in  stratified  fluids.  The  intent  of  the  in- 
verse problem  is  to  use  limited  knowledge  of  the  net  displacement  field,  such  as  might  be 
obtained  by  simple  measurements,  to  infer  characteristics  of  the  initial  disturbance.  We 
demonstrated  that  knowledge  of  the  initial  and  final  material  configurations  of  the  shape 
of  the  region  containing  the  initial  disturbance  provides  a means  for  estimating  mean 
characteristics  of  the  init>°'  kopychnal  deformation.  Thus  our  study  of  the  mathematical 
problem  suggests  that  a certain  correlation  must  exist  between  the  shape  of  the  mixed 
region  and  the  initial  disturbance. 
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